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ABSTRACT 

It  has  been  shown  that  a  finite-difference  numerical  technique  can 
be  used  to  solve  mixed  initial- and  boundary- value  problems  involving 
high-speed  elastic-plastic  flow  with  spherical  symmetry.  Numerical 
solutions  for  the  dynamic  expansion  of  a  spherical  cavity  under  a  constant 
pressure  are  presented  to  demonstrate  the  nature  and  capability  of  the 
numerical  scheme.  The  solution  for  an  elastic  material  agrees  closely 
with  the  exact  solution.  The  solution  for  an  elastic-perfectly-plastic 
material  has  confirmed  Green’s  prediction  concerning  the  motion  of  the 
elastic-plastic  boundary.  At  large  times,  the  asymptotic  solution  of 
the  dynamic  problem  is  different  from  the  quasi-static  solution.  This 
result  indicates  that  the  quasi-static  approximation  may  not  hold  in 
dynamic  plasticity.  A  non-linear  dependence  of  the  plastic  solution 
on  the  boundary  condition  is  also  observed  in  the  results. 
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1.  INTRODUCTION 


The  problem  of  high-speed  elastic-plastic  flow  with  spherical 
symmetry  is  of  interest  for  two  reasons.  First,  it  has  practical 
application  in  underground  explosions  and  in  the  detonation  of  high- 
explosives  in  solids;  secondly,  with  little  increase  in  the  mathematical 
difficulties  from  the  uniaxial  flow,  the  simple  geometry  of  this  problem 
allows  the  study  of  an  elastic-plastic  flow  involving  the  bi-axial  state 
of  stresses  and  inertia  forces.  Perhaps  this  second  characteristic  is 
equally  interesting  to  both  an  applied  mathematician  and  an  experimentalist 
who  is  concerned  with  obtaining  information  of  the  dynamic  yield  pro¬ 
perties  of  a  material  under  bi-axial  state  of  stresses.  In  the  past, 

many  authors  have  investigated  this  problem;  the  reader  may  refer  to  a 
1* 

paper  by  Hopkins  for  a  general  review  on  the  current  status  of  research 
in  this  subject.  Most  of  the  published  work  is  concerned  with  small 
deformations.  The  materials  investigated  include  the  idealized  elastic- 
perfectly-plastic  and  elastic-linear-work-hardening  materials.  However,  up 
to  present,  because  of  the  mathematical  difficulties  involved  in  treating 
two-phase  flows,  there  existed  no  analytical  solution  which  describes 
satisfactorily  all  the  phases  of  the  impact  process.  As  pointed  out  in  Section 
2  below,  only  at  the  very  beginning  of  an  impact  is  the  location  of  the 
elastic-plastic  boundary  known  and  a  simple  analytic  description  of  the 
flow  available.  However,  for  the  major  duration  of  the  impact  process, 
the  description  of  the  flow  seems  to  rely  upon  a  numerical  solution. 

The  purpose  of  this  report  is  to  show  that  this  problem  can  be  solved 
for  various  initial  and  boundary  conditions  by  using  a  finite  difference 
technique.  The  elastic-plastic  solution  obtained  for  the  expansion  of 
a  spherical  cavity  under  constant  internal  pressure  has  demonstrated 

some  interesting  facts  concerning  the  impact  phenomena.  It  has  confirmed 
1*2 

Green's  results  regarding  the  motion  of  the  elastic-plastic  boundary 
and  also  has  exhibited  the  inadequacy  of  the  quasi-static  approximation 


Superscript  numbers  denote  references  which  may  be  found  on  page  30. 
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in  dynamic  plasticity.  Numerical  solutions  of  spherical  elastic- 

3 

plastic  flow  have  also  been  obtained  by  Davids  et  al  and  by  Friedman 
et  al.^ 


2.  BASIC  EQUATIONS 


Referring  to  a  spherical  coordinate  system,  r,  0  and  <|>,  with  the 
origin  at  the  center  of  the  cavity,  only  the  equation  of  motion  in  the  r 
direction  will  not  vanish  identically  because  of  the  spherical  symmetry 
considered  here;  hence, 


do  o  -a. 

— — ^  +  2-i_i 

8r  r 


P 


a v 
at 


(1) 


In  the  equation,  r  is  the  Lagrangian  position  of  a  particle;  p  is  the 
mass  density  of  the  material  in  the  underformed  state;  v  the  radial 
velocity  of  the  particle,  and  a ^  and  a ^  are  respectively  the  normal 
stress  components  in  the  radial  and  the  circumferential  directions. 

The  stresses  are  engineering  stresses  calculated  with  reference  to  the 
undeformed  state.  This  report  considers  only  small  deformations;  hence, 
differences  between  the  engineering  and  true  values  of  stress  and  strain 
will  be  negligible. 

The  material  is  considered  to  be  elastic,  perfectly  plastic  with  a 
constant  yield  stress  in  both  simple  compression  and  tension.  Both 
von  Mises  and  Tresca  yield  conditions  reduce  to  the  form' 

a  -  a  =  +  a  (2) 

r  0-o 

The  material  will  deform  plastically  when  this  equation  is  satisfied 
and  when 

(°r  -  V  (ar  "  ae}  =  0  •  (3) 


Otherwise,  when 


or  when 


but 
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°0)  (ar  - 


<  0 


(5) 


the  material  behaves  as  an  elastic  body. 

For  the  elastic  regime,  the  material  obeys  Hooke's  law.  Expressing 
the  radial  and  circumferential  strains,  and  respectively,  in  terms 
of  the  radial  displacement  u, 

3u 

er  “  3r 


Hooke's  law  in  a  differential  form  can  be  written  as 


(6) 
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9ae  e 

8T  +  7  v  =  0 


(T) 


These  two  equations  and  the  equation  of  motion.  Equation  (l),  must  be 
solved  simultaneously  to  obtain  the  elastic  solution. 

For  most  materials,  under  moderate  pressures,  no  appreciable  plastic 
volume  change  will  occur.  Hence,  the  elastic  relation  between  the 
mean  hydrostatic  stress,  (a_^  +  2a^) /3  and  the  dilatation  au/ar  +  2u/r 


can  be  used  for  describing  the  plastic  flow.  Using  Hooke T s  law  and 
relating  to  a ^  by  the  yield  condition,  the  relation  expressed  in  a 
differential  form  is 


1  -  2v 


aa 

_ r 

at 


1/3  17  -  2/3  f  - 
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(8) 


This  equation  must  be  solved  simultaneously  with  the  yield  condition, 
Equation  (2),  and  the  equation  of  motion,  Equation  (l),  to  provide 
results  for  the  plastic  flow. 

A  convenient  representation  of  the  solutions  of  these  equations 
can  be  obtained  by  introducing  the  following  dimensionless  variables. 


x  = 


R  = 


©=  !i 


T  = 


O 

E 

cdt 


(9) 


and 


V  = 


cd° 


where  r  is  the  radius  of  the  cavity  and  c,  =  [(l  -  v)(E/p)/(l  -  v  -  2v  ) ] 
o  d 

is  the  dilatational  wave  speed  in  a*,  elastic  solid.  In  terms  of  these 
dimensionless  variables,  the  simultaneous  partial  differential  equations 
for  the  dynamic  problem  are 


2m1/2 
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and 


3R 
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©  =  R  +  1  . 

Equation  (10)  will  apply  for  elastic  deformation;  i.e.,  when 


or  when 


but 


R  -  ©  <1 


R  -  ©  =  +  1 


(R  -  ©)  (R  -  © )  <1  . 

Equation  (ll)  will  hold  for  the  plastic  flow,  when 

R  -  0  =  +  1 

and 

»  * 

(R  -  ©  )  (R  -  0)  =0  . 


(12) 


(13) 


(1M 


The  +  signs  are  for  the  situations  when  R  -  @  =  +  1  respectively. 

In  Equations  (10)  and  (ll),  only  one  parameter,  the  Poisson’s  ratio 
v  appears  in  the  coefficients.  The  equations  written  in  this  form 
provide  a  clearer  picture  of  the  dependence  of  the  solution  on  various 
parameters  of  material  properties.  For  a  stress  boundary  condition,  the 
stress  distribution  around  a  plastically  deformed  cavity  at  any  instant 
depends  only  on  the  Poisson’s  ratio  and  the  ratio  between  the  applied 
pressure  to  the  yield  stress  of  the  material. 

Equations  (10)  and  (ll)  are  two  sets  of  first  order  linear  hyperbolic 
partial  differential  equations.  Their  solutions  can  be  obtained  by 
combining  two  waves  associating  with  a  displacement  function  as  described 
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by  Hunter.  However,  this  method  requires  a  priori  a  knowledge  of  the 

position  of  the  boundary  of  the  elastic  and  plastic  regions  of  the 

1  2 

material  which  can  be  obtained  only  in  few  instances.  Green,  * 
has  proved  that  when  there  is  a  discontinuity  in  the  magnitude 

of  the  radial  strain  e  across  the  elastic-plastic  boundary,  the  boundary 

r  l/2 

will  travel  with  a  constant  velocity  c^  =  [E/3  p  (l  -  2v)]  ^  and  the 

position  of  the  boundary  is,  therefore,  known.  However,  with  the 

rapid  attentuation  of  a  diverging  spherical  wave,  this  discontinuity  in 

radial  strain  can  be  maintained  only  near  the  very  beginning  of  an 

impact.  Therefore,  other  techniques  must  be  used  to  obtain  the  solution 

which  holds  for  the  major  portion  of  the  impact. 


3,  NUMERICAL  SCHmES 

H 

In  this  report,  a  finite  difference  scheme  proposed  by  Lax  will 
be  used  to  obtain  a  solution  of  the  equations  derived  above.  The 
original  scheme  has  been  shown  suitable  for  solving  initial  value  problems 
of  hyperbolic  partial  differential  equations  in  hydrodynamics.  The 
scheme  is  extended  here  for  mixed  boundary-and  initial- value  problems 
with  regions  of  solutions  governed  by  different  sets  of  equations. 

The  results  obtained  using  this  scheme  are  compared  with  known  analytical 
solutions  for  the  case  of  elastic  waves.  The  accuracy  of  the  numerical 
results  appears  to  be  satisfactory.  An  appendix  to  this  report  describes 
the  computer  program  used  to  obtain  the  present  results. 

The  finite  difference  scheme  is  for  any  set  of  partial  differential 
equations  assuming  the  conservation  form;  i.e.. 


3F 

8x 


0  , 


(15) 


where  F,  G  and  H  are  functions  of  the  dependent  and  independent  variables. 
The  numerical  scheme  replaces  3F/9xby  a  central  difference  quotient;  i.e.. 


n 


+  1 ,k  Fn  -  l,k 


2Ax 


(16) 


12 


where  Ax  denotes  an  increment  in  x;  Fn+1  designates  the  value  of  F  at 

the  grid  point  located  at  distances  (n+l)  Ax  and  kAx  from  the  origin  in 
the  x  and  t  plane.  The  time  derivative  3G/3t  is  substituted  by  a  forward 
difference  quotient, 

Gn«.  k  +  1  ~  g^Gn+l,k  +  Gn-l,k^  ( yj ) 

At 


Finally,  the  function  H  is  replaced  by  the  average  value  at  two  neighbor¬ 
ing  grid  points  at  the  time  kA.T, 

£THn+l,k  +  Hn-l,k'  (ifi) 

The  scheme  is  equivalent  to  adding  into  the  original  equation  a  term 

2 

which  is  due  to  an  artifical  viscosity  of  the  magnitude  of  (Ax)  /2At. 

7 

The  stability  and  convergence  of  this  scheme  has  been  discussed  by  Lax. 

The  stability  criterion  is  similar  to  the  Courant-Friedrichs-Lewy 
condition,  i.e.,  At/Ax<c  where  c  is  a  quantity  which  may  depend  on  the 
variables  in  the  problem. 

The  present  equations.  Equations  (10)  and  (ll),  are  in  the  conservative 
form  and,  hence,  can  be  solved  numerically  with  this  method.  As  an 
example,  the  finite  difference  equation  for  the  first  equation  of  (10) 
is  obtained;  i.e.. 


Rn+l,k  Rn-l,k 
2Ax 


1  -  v 
2 

l-v-2v 


V  -  —  (V  + 

n,k+l  2  n+l,k 


At 


n-l,k 


Rn+l,k  ® n+l , k  Rn-l,k  ®n-l,k  _ 

*  >  U  ■ 

xn+l,k  xn-l,k 

After  being  converted  to  finite  difference  equations,  the  three  equations 
in  (10)  and  in  (ll)  will  provide  three  algebraic  equations  sufficient 
for  solving  the  values  of  the  dependent  variables  at  t  =  (k  +1)At; 
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from  the  known  data  at  a  previous 


namely,  R,(k+1,  Vn>k+1  and  @,jktl 
time  kit;  namely,  Rn_1>k,  9n+lk,  ®n_1>k,  V1>k  and  At  each 

point,  the  finite  difference  equations  for  elastic  deformations.  Equation 

(10),  will  he  first  used  to  obtain  R  ,  ..  and  0  .  .  .  If  the  difference 

n,k+l  n,k+l 

is  greater  than  unity,  then  the  calculations 


between  them  R  ,  , _  -  0  ,  . , 

!  n,k+l  n,k*l 

will  be  revised  by  using  the  finite  difference  equations  for  plastic 
flow  which  are  derived  from  Equations  (ll). 


The  above  scheme  can  be  applied  very  conveniently  at  points  other 
than  those  on  the  boundary.  However,  when  the  central  point  is  on 
the  cavity  surface,  the  point  xn  ^  ^  will  be  outside  the  body,  and 
fictitious  values  must  be  assigned  for  the  dependent  variables  at  this 
point  in  order  to  use  the  Lax  numerical  schemes.  The  present  investi¬ 
gation  shows  that,  for  this  situation,  satisfactory  results  can  be 
obtained  by  using  a  different  numerical  scheme.  In  the  new  scheme,  the 
Lax  method  is  applied  in  reverse  to  those  dependent  variables  whose 
values  are  not  specified  on  the  boundary.  Accordingly,  a  central  differ¬ 
ence  formula  is  used  for  the  time  derivative  and  a  forward  difference 
formula  for  the  space  derivative.  As  an  example,  the  space  and  time 
derivatives  for  the  dependent  variable  V  are 

V  -  V 

3V  =  n,k+l  n ,k+l 

8T  '  24t  (19) 

V  -  —  (V  +  V  ) 

av  =  n+l ,k  2  n,k+l  n,k-Iy  . 

3x  Ax 


However,  for  the  derivatives  of  the  dependent  variables  whose  values  are 
specified  on  the  boundary,  R  in  this  case,  the  Lax  scheme  is  kept.  As 
a  result,  a  fictitious  value  ^  must  be  assumed,  which,  together  with 

the  other  two  unknowns  V  .  _  and  ©  .  ,  can  be  obtained  from  the  three 

n,k+l  n,k+l 

finite-difference  equations  derived  from  Equations  (10)  or  Equation  (ll) 


using  Equation  (19).  One  of  the  unknowns  in  the  calculation  at  an  interior 
point,  jr+1  is  now  a  known  quantity  which  has  the  specified  value 

of  R  on  the  boundary  at  x  =  (k  +  l)Ax. 

k,  RESULTS  FOR  ELASTIC  CASE  -  ACCURACY  OF 
THE  NUMERICAL  SCHEME 

The  case  of  a  constant  pressure  applied  suddenly  on  the  surface  of 

a  spherical  cavity  within  an  infinite  elastic  solid  is  studied  to 

demonstrate  the  capability  of  the  numerical  scheme.  In  Figures  la  and 

lb,  the  principal  shear  stress  in  the  present  numerical  results  is 

8 

compared  to  those  obtained  analytically  by  Hunter  and  numerically  by 
o 

Chou.  Chou  uses  the  method  of  characteristics  in  obtaining  his  solution. 
Considering  the  rather  large  grid  size  (twenty-five  points  per  cavity 
radius)  used  to  obtain  the  present  solution,  the  agreement  between  this 
solution  and  the  other  solutions  is  satisfactory. 

However,  near  the  elastic  wave  front  where  a  jump  in  the  values  of 
the  shear  stress  exists,  a  large  discrepancy  is  observed  between  this 
solution  and  the  exact  one.  The  discrepancy  is  mainly  due  to  the 
artificial  viscosity  which  is  associated  with  the  numerical  scheme  and 
tends  to  smooth  the  sharp  discontinuities  of  the  solution.  Since  the 
artificial  viscosity  is  equal  to  (l/2)  (Ax)  /At,  the  effects  can  be  greatly 

reduced  by  increasing  the  ratio  of  Ax/ Ax,  as  indicated  by  comparing  the 
results  obtained  with  At/Ax  =  1  and  1/2.  For  this  case,  Ax/At  =  1  is 
the  maximiam  value  allowed  by  the  stability  criterion;  at  this  Ax/At 
ratio  some  small  oscillations  in  the  solution  which  are  not  shown  in 
the  figures  have  already  occurred  near  the  wave  front.  The  effects  of 
the  viscosity  can  also  be  reduced  by  diminishing  the  grid  size  as  seen 
in  the  results  for  Ax  =  0.01  and  0.02.  However,  considering  the  increase 
of  computing  time,  the  second  method  of  reducing  viscosity  is  not  as 
efficient  as  the  first  one.  Roberts10  has  proposed  a  much  more  effective 
method  for  improving  the  accuracy  of  this  numerical  scheme.  The 
technique  involves  combining  several  solutions  of  first  order  accuracy 
obtained  with  a  fixed  Ax/At  ratio.  Solutions  with  accuracy  of  order 

greater  than  one  can  be  obtained  very  easily.  IBCOTRTY  0?  ILS.  AR NTi 

SIT TTF0 
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Dimensionless  Maximum  Shear  Stress,  (R-8)/R 
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Figure  lb.  Comparison  of  numerical  and  exact  solutions  tor  the  elastic  expansion 
of  a  spherical  cavity;  the  variation  of  maximum  shear  stress  at  the 
cavity  surface  after  the  sudden  application  a  constant  pressure, 

R  =  1 ,  on  the  surface, 
o 
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5*  RESULTS  FOR  ELASTIC-PERFECTLY-PLASTIC  MATERIAL 


The  numerical  scheme  just  discussed  has  also  been  used  to  obtain 
the  elastic-perfectly-plastic  solution  for  a  spherical  cavity  expanded 
under  a  constant  internal  pressure  which  is  suddenly  applied  at  t  =  0. 
Typical  sets  of  results  are  shown  in  Figures  2a  through  3b  for  the 
cases  in  which  the  magnitudes  of  the  internal  pressures  are  2  and  20 
times  of  the  constant  yield  stress  of  the  material.  In  the  results  for 
the  latter  case,  the  evidence  of  two  discontinuities  can  be  seen  in 
the  profiles  of  stresses  along  the  radial  axis.  The  one  which  is  farther 
away  from  the  cavity  surface  is  the  elastic  wave  front  while  the  one 
propagating  behind  is  the  plastic  wave  front  which  separates  the  regions 
of  plastic  and  elastic  deformation.  The  constant  speed  of  propagation 
of  the  elastic  wave  front  agrees  closely  with  the  exact  value  of  the 
dilatation  wave  speed,  c^  obtained  here.  The  magnitude  of  the  jump  in 
the  dependent  variables  across  the  elastic  wave  front  also  seems  to 
agree  with  the  exact  solution  which  is  represented  by  the  dashed  curves 
in  Figures  2a  to  2c. 

The  elastic-plastic  boundary  moves  at  a  constant  speed  when  a  dis¬ 
continuity  in  the  radial  stress  across  the  boundary  is  visible  near  the 

beginning  of  the  impact.  The  value  of  the  speed  agrees  with  the 

1  1/2 

theoretical  one  predicted  by  Green,  namely,  c  =  [E/3p(l  -  2v)] 

1? 

This  result  is  shown  in  Figure  b  in  which  the  position  of  the  elastic- 
plastic  boundary  at  all  times  is  plotted.  Figure  b  also  indicates  that, 
for  the  present  boundary  condition,  the  movement  of  the  elastic-plastic 
boundary  will  slow  down  as  the  discontinuity  across  the  boundary 
disappears  because  of  the  geometrical  dispersion.  As  long  as  the  pressure 
on  the  cavity  surface  remains  unchanged,  the  motion  of  the  elastic- 
plastic  boundary  will  come  to  a  halt  and  then  reverse  its  direction 
towards  the  cavity  surface.  This  result  is  quite  different  from  that 
observed  in  the  static  solution  for  the  same  elastic-plastic  problem; 
when  the  constant  pressure  is  applied  quasi-statically  at  the  boundary, 
a  plastic  zone  next  to  the  cavity  surface  is  always  maintained  by  the 
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Figure  2a.  Numerical  results  for  the  elastic-plastic  expansion  of  a  spherical 

cavity;  the  distributions  of  radial  stress  along  the  radial  axis  for 
various  times  after  the  sudden  application  of  a  constant  pressure, 

Ro  -  2,  on  the  cavity  surface.  Arrows  indicate  the  positions  of  the 
elastic-plastic  boundary. 


Dimensionless  Circumferential  Stress 


Figure  2c.  Numerical  results  for  the  elastic-plastic  expansion  of  a  spherical 

cavity;  the  distributions  of  particle  velocity  along  the  radial  axis 
for  various  times  after  the  sudden  application  of  a  constant  pressure, 
R0  =  2,  on  the  cavity  surface.  Arrows  indicate  the  positions  of  the 
elastic-plastic  boundary. 
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Figure  2d. 


Numerical  results  for  the  elastic-plastic  expansion  of  a  spberleal 
cavity;  the  distributions  of  radial  displacement  along  the  radial  axis 
for  various  times  after  the  sudden  application  of  a  constant  pressure, 
Ro  =  25  on  the  cavity  surface.  Arrows  indicate  the  positions  ot  tne 

elastic-plastic  boundary. 
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Figure  3a.  Numerical  results  for  the  elastic-plastic  expansion  of  a  spherical 

cavity;  the  distributions  of  radial  stress  along  the  radial  axis  for 
several  times  after  the  sudden  application  of  a  constant  pressure, 

Ro  =  20,  on  the  cavity  surface.  Mrrows  indicate  the  positions  of  the 
elastic-plastic  boundary. 
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Figure  3b.  Numerical  results  for  the  elastic-plastic  expansion  of  a  spherical 
cavity;  the  distributions  of  circumferential  stress  along  the  radial 
axis  for  several  times  after  the  sudden  application  of  a  constant 
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applied  pressure.^  Differences  "between  the  static  solutions  and  the 
asymptotic  dynamic  solution  at  large  times  can  also  he  observed  in  the 
distributions  of  the  stresses  and  displacements  along  the  radial  axis 
(Figures  2a  to  2d).  A  larger  deformation  and  a  higher  radial  stress 
occur  in  the  dynamic  case  than  do  in  the  static  case.  These 
differences  are  not  likely  to  be  attributable  to  errors  in  the  numerical 
solution,  because  the  comparison  of  the  results  obtained  with  two 
grid  sizes  does  not  suggest  the  possible  existence  of  such  large  differences 
in  the  solution  (Figure  2a) .  The  maximum  radius  of  the  elastic-plastic 
boundary  is  also  larger  in  the  dynamic  case  than  in  the  static  case. 

These  disagreements  between  the  dynamic  solution  at  a  large  time  and 
the  static  solution  are  probably  due  to  the  residual  stresses  produced 
by  the  inertial  force  occurring  in  the  dynamic  deformation  process. 

This  result  implies  that,  in  some  circumstances,  the  quasi-static  approxi¬ 
mation  may  not  be  valid  in  dynamic  plasticity.  The  quasi-static 
approximation  in  which  the  static  solution  of  a  problem  is  used  to 

describe  the  long-time  behavior  of  the  dynamic  solution  holds  in  many 

8  11  12 

situations  in  dynamic  elasticity.  ’  9 

Figures  h  and  5  show,  for  all  times,  the  dependence  of  the  particle 
velocity  and  of  the* position  of  the  elastic-plastic  boundary  on  the 
magnitude  of  the  applied  pressure.  Although  all  of  the  results  for 
various  applied  pressure  seem  to  be  similar  qualitatively,  the  magnitudes 
of  various  parameters  in  the  solution  seems  to  be  quite  sensitive  to  a 
change  in  the  pressure.  The  dependence  is  rather  nonlinear.  It  is 
also  interesting  to  note  that  the  nonlinearity  also  appears  in  the 
result  of  the  total  duration  of  the  dynamic  deformation,  namely,  the 
duration  during  which  the  acceleration  of  the  material  is  not  vanish  ingly 
small.  This  duration  for  the  elastic-plastic  problem  increases  with 
increasing  applied  pressure,  while,  in  the  elastic  case,  the  duration 
is  a  constant  as  shown  in  Figure  5*  Accordingly,  the  initiation  of 
slight  plastic  deformations  tend  to  change  considerably  the  nature  of 
the  flow. 
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6.  CONCLUSIONS 


It  has  been  shown  that  the  governing  equations  for  an  elastic- 
perfectly-plastic  flow  with  spherical  symmetry  can  be  solved  through 
the  use  of  a  finite  difference  technique  whose  nature  is  rather  well 
understood  in  comparison  to  other  similar  methods.  With  a  slight 
modification  of  the  original  scheme  proposed  by  Lax,  the  method  can  be 
applied  to  solve  a  wide  variety  of  mixed  initial -and  boundary -value 
problems  involving  elastic  and  plastic  flows.  The  comparison  of  the 
results  for  the  dynamic  deformation  of  an  elastic  cavity  obtained  using 
this  method  and  an  exact  analytical  method  shows  that  the  accuracy  of 
the  numerical  solution  is  acceptable,  except  probably  in  the  region 
where  a  discontinuity  in  the  dependent  variables  appears.  Large  errors 
may  occur  in  this  region  because  of  the  artificial  viscosity  added 
into  the  equations  by  using  the  numerical  scheme.  Accordingly,  for 
solutions  containing  discontinuities,  an  investigation  of  the  error 
is  necessary  by  varying  the  grid  size  and  Ax/At  ratio.  If  high  accuracy 
is  required  for  the  solution  near  the  discontinuities,  it  is  advisable 
to  locate  and  obtain  the  discontinuities  in  the  solution  using  an 

U 

analytical  method  as  shown  by  Friedman  et  al. 


For  the  dynamic  expansion  of  a  spherical  cavity  in  an  elastic- 

perfectly-plastic  solid,  the  numerical  results  for  a  suddenly  applied 

internal  pressure  confirm  Green* s  analytical  prediction  that  an  elastic- 

*  1/2 

plastic  boundary  propagates  with  a  constant  speed,  c  =  [E/3p(l  -  2v)]  . 

P 

when  discontinuities  in  stress  and  particle  velocity  occur  across  the 
boundary.  Within  the  major  duration  of  an  impact  the  elastic-plastic 
boundary  moves  with  a  variable  speed,  of  which  the  magnitude  cannot  be 
predicted  analytically.  Therefore,  a  numerical  technique  seems  to 
be  unavoidable  in  the  investigation  of  spherical  elastic-plastic  flow. 


The  same  results  have  also  shown  that  the  quasi-static  approximation 
which  is  often  used  in  dynamic  elasticity  may  not  necessarily  hold  in 
dynamic  plasticity.  The  solution  for  the  dynamic  problem  investigated 
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does  not  approach  the  static  solution  as  a  limit  when  time  elapses. 
This  result  is  probably  due  to  the  residual  stress  produced  by  the 
inertia  force  of  the  impact. 

As  far  as  magnitudes  are  concerned,  the  results  show  a  rather 
nonlinear  dependence  of  the  elastic-plastic  solution  of  the  internal 
pressure.  This  nonlinearity  seems  to  indicate  that  the  nature  of  the 
flow  would  be  changed  considerably  by  the  initiation  of  a  slight 
amount  of  plastic  deformation. 
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APPENDIX 


A  COMPUTER  PROGRAM  FOR  THE  DYNAMIC  EXPANSION  OF 
A  SPHERICAL  CAVITY  BY  AN  INTERNAL  PRESSURE 
(ELASTIC-PERFECTLY-PLASTIC  SOLUTION) 

The  computer  program  used  for  obtaining  the  present  results  is 
attached  at  the  end  of  this  Appendix.  The  program  was  written  in  the 
FORAST  language  by  Mr.  T.  Addison  of  the  Computing  Laboratory.  (The 
reader  may  consult  one  of  the  references  at  the  end  of  the  Appendix  for 
a  description  of  the  computer  language. )  The  program  can  be  executed 
with  the  BRLESC  Computer  which  is  currently  being  used  in  the  Ballistic 
Research  Laboratories. 

The  computer  program  is  based  on  the  finite  difference  schemes 
described  in  the  main  text  of  this  report  and  will  be  able  to  provide 
information  of  stress,  displacement,  and  velocity  around  a  spherical 
cavity  which  is  being  expanded  by  an  internal  pressure  specified  with 
card  numbers  157  and  171 •  The  pressure  at  the  boundary  is  shown  as  a 
constant,  RR,  but  can  be  changed  to  a  function  of  time,  f(x);  t  is 
represented  by  T  in  the  program.  The  physical  quantities  R,  V,  and  U 
are  represented  by  R,  Q,  V  and  U,  respectively. 

The  program  can  be  used  for  problems  with  non-vanishing  initial 
values,  if  the  section  between  cards  13  and  39  is  properly  modified. 

This  section  is  now  written  to  represent  the  similarity  solution  near 
t  =  0  for  this  problem.  The  similarity  solution  was  put  in  the  program 
as  an  alternate  way  of  calculating  the  results  but  was  omitted  in 
obtaining  the  final  results  presented  in  this  report. 

The  numerical  calculation  proceeds  in  increments  of  time,  At,  while 
at  a  given  time,  t,  it  begins  with  the  two  data  points  nearest  to  the 
boundary,  x  =  1  and  shifts  in  the  positive  x  direction.  Calculations 

a  +  -nn-irvf-Q  fhan  f.ho  hmmHay»v  -nr^i  rvh  ahou-n  in  r»  q  c:  +•.  r\  70  •  ■i  n 

cards  60  to  65,  the  plastic  equations  are  used,  and  in  cards  72  to  79, 
the  elastic  equations  are  applied.  Computations  at  the  boundary  point 
occur  only  once  every  other  time  step,  as  shown  in  cards  155  to  170  for 
an  elastic  case  and  in  cards  171  to  182  for  a  plastic  case. 
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The  input  data  concerning  material  properties,  specifications  for 
outputs  and  the  dimensions  of  the  grid  used  for  calculations  are  punched 
in  the  two  cards  at  the  end  of  the  program  which  will  be  read  by  the 
machine  upon  the  order  written  in  cards  9  and  11.  The  input  data  include 
the  following  information: 


Program  Notations 


vs 


Physical  Quanitites 


DX; 

DT: 

DTT: 

DTI: 

TMAX: 

RR: 

NU 

RTO: 

A: 

TO: 


Ax 

At 

t  of  the  first  output 

interval  of  t  for  output 

t  of  terminating  calcu¬ 
lation 

Radial  stress  at  the 
boundary,  x  =  1 

v 

initial  value  of  radial 
stress 

a 

initial  time. 


The  output  of  the  calculation  will  be  stored  in  a  magnetic  tape  and 
can  be  tabulated  in  a  numerical  form  (cards  112  to  126)  as  well  as  plotted 
in  a  graphical  form  (cards  to  lUT ) -  The  format  of  printing  is 
indicated  in  cards  l6  and  IT.  In  plotting,  the  distributions  of  stresses, 
particle  velocity  and  displacement  along  the  radial  axis  will  be  presented 
by  color  curves  in  a  single  graph  for  sach  time,  t.  The  size  and  scaling  of 
the  graph  is  controlled  with  cards  No.  128  and  129*  The  curves  of  R, 

@,  V  and  U  are  all  referred  to  the  same  scale  which  is  chosen  for  R,  and 
hence,  it  is  sometimes  necessary  to  adjust  the  values  of  V  and  U  with 
proper  numerical  factors  before  plotting,  as  shown  in  cards  132  and  1^9* 


32 


'  .  1 

L 

T - HTXTT751  H1N - T 

X  MAXOl930UU)LlNt.S  4 

BLUCir-Y30UO)AHl-AH444  )UO-U5oUO  )  V u - V5 0 0 0  > R 0 -R50 0 0 ) 0 0 -Q5 U 0 0  )  5 

tiL0Cill*il2)Dl-U2>  b 

bYN  (Y=Xl>  '  ‘  / 

-TART  NEAD-(!-W2)-lUl)D2)  b 

wEad(ux)dt)uit)dii)tmax)rh)  v 

t.N  i  E«  l  sEr  .  I  0  )  8  )  0  )  U  )  100  >  10 

"RXA'D  rWj)  xfTH  1 0  )  A '  )T(j  )  "  .  il 

«EAD-(Fl)-(12)N0b.AT(U)X  12 

'TST  FFJ I  < «.  =  i )  *  fNT(KK=l2*l)X  X  =  X1*  lNT(N=4j*  t  =  T  U  X  14 

U2U2=U-NU-2*NU**2)/(1-NU>X  R=R10X  CUCU=buH  i  (  C2U2  i  14 

0N=<i-i\iU)/4ll*Nu>(l-2*NU)**2X  CN=bQK I  (CN)  15 

IFlRU-NU/tl-NU)  )>  1)GOTO<3.0>  lb 

1.0  «1»I=KI0X  U1  *  I  =  <N|j/(l-Nu)  )RX  vl»  I=1’A'*K*C2U2  I/' 

mm=ia  lb 

COuNTT -3)  lM(  HGUTQ ( 1  ,“5Tx  iNT(F?Il  iV 

uOiO(HjN)  cu 

T7u  iF  ( xJU  m*to  >  liu  i  U(  4 .  (i )  cl 

Hi,  1  =H  1  OX  Ul,  1=R1,  I -IX  _  _  _  _  4  2 

V 1 »  Is  -  A  •  (  ( 1  +  N  U  )  -  <  >5  ( 1  -  N  U  )  *  *  37(1  +  MU  )  (  i  ”  2  *  N  U  )  *  *2  )  **  .  5  e  4 

CUNT  ♦(  3^  1-NU)  /  ( 1+NU)  )  «*  .5«Rl  >  I  _ _ _  44 

X a X  +  2 *  U X X  MM  =  2 . 0 X  COUN  T (0  > I  N (  V )  UU  I  0  i 3.0)  .  £5 

OOlO(HuN)  4b 

A  .  0  Hi*  Is(l"NU)/(l_2*NU)X  Ul,  l=M/(l-2*NU)/. / 1 ,  i  =  -  (  1  ♦  h  u  J  a  '  e.  / 

XaX  +  2*UXX  MM  =  3.0X  £t! 

l-OUNTT  oJ  I  nTTT  g  u  T  0 "(  3 . 0  )"  .  ”  "  “  £V 

HUN  CLLAKt 5)N0b. A ! (R1 >%  CLt AR ( 5  )  M-S . A  I  ( V ±  )  X  CLt AK ( 5 ) Nib . A T ( Ul )  ou 

huN  FUNcflcnri t I al  condYt ionS>x  ex  i eh (Hunch  u )  '  '  "  ”  01 

HUNCH<  J  X  V  K  >  U2 

cwyt  c  T> 

UEI  ( 1  =  0) J  =  1)X  CC1  =  D1X  CC?=D2  u 4 

HWl  "  PuNCH-(Fw5>-rCCi  )CC2)  J)X)V1,  I  IR1,  1)01,  i)j)  05 

uCl  =  CC2  =  0  s  fc 

xixTi •Uxx  Tnc (  jI'jTjT) x“couKWriN ( i')'iiut dT y uni V  0/ 

_ bE  *  <  1  =0  ) _ 0B_ 

HUNCH<tJEC>  Pa^E>  ov 

COHH _ EO^PuTE  CONSJaJnTS  4U 

C  C  N  (j  s  (  1=NU_2*NU**2  )  41 

Arl/(l-NU-«i*NU*NU  >  A  *  X  2DT  =  l/2*Dr%  lN^l-NU  42 

Si 1 a2U| *iN  *  A  X  bl 2  =  1 /2 • U  X- 1  %  E£ 1=  NU/DX"iX  b23  =  i  / A*uT  +0 

_ B31«A(  ■  5*1N/UX-NU  )X  832» -2  U  T _ _ _ J*4_ 

A 1 1  a2l)  I  X  A22=«20T%  2Dx«1./2*EX%  A12  =  A  •  ,C2U2*2Dx  '45 

A2l»(2uX~l)  ( 1  /  3  *  A '  (1-2»NIJ))  _  __  _ _ _  _  __  ‘*6 

?E]“(BY  =  l”)  ‘  ’  ’  ‘  . .  ’  . .  '  "  '  ?/ 

4.1  _ J^1  =  1*_N2  =  2X  l=_ltur  X.S.b  T  (  1  =  0  > _ _ ...  .  4b 

1NC(K  =  X*D  4  V 

_ (F-1  Nl  (Ks2(K/2)  )  GQTQ(tVEN)X _ _ _  __~±u 

V K 2 a  V  1  X  QK2  =  Ul%  HKpaRl’  X  SET  (  A  Yal  )  GO  I  0<  f  •  2  )  5l' 

EVEN _ tN  I  (  N  =  N  ♦  j  )X  b  F  I  (  A  Y  =  2  ) _ _ _ _ _  __  _  ?2 

Y K 1  a  V 1 X  QK 1 s U 1 X  H K 1  a  R1 %  yK 1  a  C 1  ”  .  04 

4.2  U0l0*AY(AYl)AV2)  54 

A  Yl  XR  =  x  1  ♦  il2 * D XX  N 2  =  N 2  +  2 X  O0 T E  ("4 .5  )  .  55 

A  Y 2  XN  =  X1*  N1»PXX  Nl  -Nl  +2 _ =b_ 

A .  5  C1  =  D  I((1»2*NU)(V2*!-V1»I)/2«I>+1nu-.5)iV2»1  +  vi,1)/XN)/  5/ 

CONT  CCNU*A,  +  .5(H2*  I  +  Ri ,  I  ) - , 5 ( 02  *  l  +  Ul,  1  )  5b 

1  F-AbS ( C 1< 1 .  ) UO 1 0 ( EQ2  )  . .  '  5V 


33 


.  IUHH  tQQAIlON  1  PLASTIC  .  _ 

lv  *DI  •  A ' *C 2 D 2  <  (H?,  I-Rl,  1  )/2*DX*’2/XN)*(v2»  l  ♦  V 1  >  i  )  .5 
Kl.»  I*DI  (  (Viti.  i-Vl,  I  )/2*DX+<V2-  1  *  V 1 ,  I  ) /XN  ) /i  ( i-2*Nu  i  A  '  ♦ 

CONt  ,B(K2.I*R1 ,  1  ) 

_  Ul>  l8*!,  1-1 _  _ _ 

01 » I*<0T(V2,  I*V1,  I )*U2. I+Ll,  I  > .5 

_ _ yi_Lis,y_ _ 

l  P2  (  U!?U  )  <  O^L  )  (VI V  I  )  XX  0  GOEi  TO  BII  y  UF  Vl 

_ OQMH  ZtRO  GOES  IU  BIT  9  Of  V1=PLABT1C  EQ  ■ _ 

COUN  t  <  N ) I N l 1 >GUTO (4.2) 

_ tr-lNHK«2U/2>)(»0T0(B.8)X  GC  1  0  (BOND  ) 

COMM  EQUATION  2  ELASTIC 

LQ2 _ iy  =  C2L)2«A,»IJI((R2.I-Rl,I)/2»LX*(B2.  l*Rl.  1-02,  I-Ul  ,  1  ) 

CONT  /Xfy))  ♦  (  V2»  1  *  VI  ►  I  )  .  5 

K1 »  ISL)(  ( IN  (  v  2  ,  1-V1,  I  )/2*DX+NL(V2»  l  +  V 1 .  I  )  /  XN  )  /  (  CONC  )  A  - 
CONT  ♦.5<H2,  I+Kl,  1  ) 

_ ui»  l*o  i  (  n  u  ( v  2 ,  t « v  1 ,  i  )/2»nx*(v«i.i*vi»t)/2«XN)/  c  c n  u»  a  ; _ 

CONT  ♦ . 5 l U2  » I +U1 1  1  ) 

_ 0 1  >  l  =  (liT  (  V  4  ,  1*»1»  1  )*U2>  I  ♦  C  1 .  I  )  .? _ _ _ _  _ 

VI.  I  *  I  v 

_ IP2(V4l) (Q4l)  (Vj,  1 )X _ 

comm  one  goes  10  bit  9  or  visEsasiic  to. 

_ _ COUNI  (N)  IN(  I  JQUT0(4.2)% _ _ _ _ 

iF-IN1 (K»2(K/2) )G0T0<8.8>*  GC  IO(0ONU) 

B.8  j>EI  (  1SU)  J=1)X  XsXlX  CC1»D1X  CC2=D2 _ 

i  F  -  I  N  r  (  K  >  H)  GO  T  0  (  9 . 5  ) 

_ tr-INT (K  =  2TK/2) ) got  0  <  4 . 1 ) _ 

I^UNCM- (  FW1 ) -<■  T1ME=  >T<  CYCLE  NO.  >(K  j 

_ t-N  I  E*  (  r>UNCM  B  ) _ 

rUNCH<  I  X  V  H  > 

_ CONTC  "  0 _ U _ [> _ _ 

9.0  »HX(Vl, I ) ( «0 ) (MM) 

_ COMM  BIT  NINE  GOtS  rO  RIT  ONE _ 

B8 ( MM ) ( 04L  )  ( MM ) 

_ comm  bit  one  QNcr  results  in  mm«o  oh  i _ _ 

I  P2  (  ubu  )  ( (HL  )  (  VI ,  i  )XX  REPLACES  d*  T  9 

_ runch-(Fh5)-(cci)Cc:2)  j)x> vi«  PRi.noi, 1)01,1) mm)  _ 

X=X*2*UXX  CC1*CC2=0  X  1  NC  T  J"v  +  1 ) 

_ LQUiyi  (y*i)  1  N  (  I  )  G  U  1  0  (  9 , 0  )  X  PUncekEoEC  '  kAGE> _ 

co 1014.1) 


oU 

01 

C 

bb 
0  4 
bb 
Ob 
0  7 
ob 
ov 
f u 
/ 1 
/'£ 
/  b 

/  H 
/  6 
/ t 

t  / 

/  b 

/  v 

QU 
01 
0  2 
ob 
04 
00 
ob 
0  / 
ob^ 
0  v 
vu 
vi 

1A. 

y  * 
y  4 

vo 

v  b 
v  / 
y  0 
y  v 


v  ,0 


AT-  i  INI -NO  TCK  =  kK)  GO  to  ( 10,1  ) 


AF-N^l  l  T  =  U  I  I  )  W  I  T  H  t  m ( d  f 
UT  I  »U  I  1  >dt  1 _ 


)G0TC^1U.1) 


iF-INI-NOTtK=2(K/2))GOTO(HuNl 
iClsU^spX  X  =  X  ♦DX _ 


luU 
lul 
I  u  2 
1  U  b 
lu4 


ryNCH- ( FW1  > 
t-N  I  E*  CRUNCH 


’<  I  I  ME  = 
B) 


> T <  CYCLE  NO.  ><<  ) 


KUNCH< 

cqnt ^  ~  J _ y_ 

OHX ( VI,  I  )  ( K8  >< MM ) 


8.V 


iu5 

lub 
i  u  / 

~iuy 

iiu 


iP2(U!»0)(04L)(Vlfl) 

rUNCH-(FW5)-(CCl)Cc2) J)X)V1, I)Rl,I)Ql.i)Ul 

.  1  )  MM  ) 

111 

112 

X  =  X*2*L)XX  i  NC  (  J  =  J  +  -]  ) 
c 0 U N I  IN  ) 1N( I )GOTo<8.9)\  GOTl(PER) 

lib 

li4 

A  N  (  t  kK  =  KK  ♦  l  2  > 

110 

HgNl 

RlJNCM-(Fwl)-<-  1  I  ME  =  >  T  <  C  *  CLt  NO.  >(R  JX 

LNTEK(EUNCH  B) 

lib 

RUNCM<  I  X 

CONT  ^  U  U  I> 

V  H 

->  11/ 
lib 

V  .  / 

9HX(V1,  I ) (KB)  (MM) 

1 1 V 

_ 8.1  T.-ftilNE  B.Ufct>  Tit. B1I— QN£ _ 1A1L 

bS(MM) <04L> <MN>  141 

COMM  BIT  UVE  (INLY  RESULTS  IN  MM*0  OH  1  1C4 

1  P4( U5Q  )  (  0  AL )  (VI.  1  )  XX  REPLACE*  BIT  *  TZT 

_ HUNCH»(FW5)-(CCl)Cr2)  J)X)V1>  P  R1 , 1  )  01 ,  1  >  ul ,  1  >  MM  ) _ 144 

X«X*4*UXX  CC1»CC2=0  X  1NCU»V  +  1)  145 

_ COUN  I  (N*l)  1N(  1  )G0  I  0(9.  /  )X _ 14C 

WF  P(JNCH<ljeCI  HaGE>  14/ 

_ XMAXg14.QX  HMAXsl^.X  XMIns1.QA  HMlN=“4.3 _ ICC 

HX=RS  =  KMAX/ 10X  XX  =  xS=(XMAX-XMN)/46  14V 

6.6Q _ *E  I  ( I =U  )X _ 10 U_ 

Y  *  X 1  101 

6.61  _ Vj.  I=-vl,  1/5.QX  Y1  ,  I=Y>  I*Z»DXX  01.  13-Ui>  1/3UQ. _  10  4 _ 

COUN I (N*l  )  INI  I  )GOTO<6.61 )  103 

_ C0I0>BY(AY*)(AY6)(AY3)(Ay4) _ 104 

A  Y5  *E I (XB=4) ( YB=14 )X  133 

AY/ _ 1NC(BY=BY*1)0QIQ(9.8) _ _ _ _ _ _ 10C_ 

AY6  *E  I ( XB  =  3l ) 00 1 0 ( A  Y / )  ’  10/ 

A  YO _ *E  I  (  YH  =  42  )  CU  I  Q  (  A  Y  /  ) _ _ _  _  IOC 

A  Y  4  *E I  ( XB  =  4 )  ( b  Y  =  1 ) X  10  V 

V.b  bTUR _ 1  a  (j 

HUNCH<XB=  >Xb  <  YB=  >  YB  1  a  1 

_ _  tNIEH(HLOT  •  *  >  1 )  >Xb)  >YB)05l'05U)XS)RS)AHl)  )  _  __  ___  1A4 

t-N  I  EH  (  -‘LOT  .0)1  )1  >1  )3>  Y  >Ql  )  ,  1  )  1A0 

_ _ tN  I  EHt^LOT  .0)1  )1  )1  )4)  Y)Rl  >  .  1  )  _  _  _ _  1A4 

tNIEHlHLot  .0)1)1) 1 ) 5  >  Y ) U 1 > .  I  )  1  A3 

_ tNIEH(HL0T.0)l)l)l)6)Y)Vi  ),  1) _ ’ _ 1 A t 

LNIEH(-'L0T.A)1)XX)RX)XMIN)XMAX)HM1N)HMAXJ  1A/ 

_ *E  I  (  i  =  U  ) _ _  _ _ _  lAb 

6.6 2  Vl, l=-vl, 1*».0X  Ui,  1=-Ul,  i*50o  1  a v 

_  CQUN I  ( N*l )  1N( I )G0 1 Q( 6.62)  J _ _  _  130 

V  .  V  1 F -  1 N I  ( B  Y  > 1  >  U  0  I  0 ( 1 n . 1) X  E  A  T E  fi l P  L  6  I  . K > 1 1  131 

IQ  .  1 _ IF  (  T  >  I  M  Ax)COIU(NHkqB) _ 134 

1F-1NI  ('KsIOHiOIO(NPROB)X  C0TC(4.l)  130 

NBKUB _ LN  I  E*  (  lNDPLU  I  )  )  l )  X  GOTO ( N . PRCc ) _ _  _  _  _ 13  4 

bUND  *E I ( 1 =U )  133 

_ _ blA«4Ur»iN«A»vKl  +  t  . 5 / D  X  » 1  )  RK2-QK1  _  13^ 

?Q  =  R*  13/ 

_ b24sUKl/A*U  I  *P*NU*vK2/UX-(  NU/LX-1 )  VKl _ 1 3_e_ 

b3AB-Hu/DT  +  4UT*RK4  +  A<lN*vK2/nx-(.3*lN/ijX-Nu)VKl)  13V 

_ YQ  8(bO2(B40»Rl4-Bg4)-B34«B23*B12) _ / ( ti62 C b23*Bl 1 -B21 >  _  _  loo 

CONT  -b31*d23*bl4)  “  id 

_ UO _ -  (  a24-b41«vc] )  /r23 _ _ _ _  104 

i F  ~  A b b ( R o - U 0  > 1  ) 0 u  T  O (BO NO i  )  103 

_ _ U0  =  (  VU  +  VKl  )  U  I  »UKl _ 1  o  A _ 

MOvElN<NOS.fKOM<VO,N/-l)TL  (vl.N/”l)  105 

MOVE(N)NOS.l-wijM(RO,N/-l)TC.JRl,N/-i)  _  _  _  ice 

MOVE  IN)  NOS  •  ►  NOM<  UO  ,  N/-1 )  TL  (C1.N/-1)  "  ~  10/ 

_  _  nQV£(N)NOS.K*l)M(UQ,N/»l)TL(lll.N/-l)  _ _  _  _  __  lot 

'  "v'i.=Yux  oi  =  uox  Ri  =  Rox  ui»uu  .  "  10V 

_ up  i  o  ( a .  e ) _ _ i  t  u_ 

BONDI  HO  =  RH  '  “  l/'l 

_  _ A20°;t*Q-.3*NK4)/UT*(VK2/i;X»UDX-l)V*l)/L>»A’(i-4«NU)  _  __  1/4 

Alb=40I *vK1*A' «C4U? ( 2DX*RK2*2 >  ’  1/0 

_  ^0 _ f  <  a22*aiowai4*a232/ <  All^Ac2-A41*Aij/ )  _  _ _  l/A 

.  U  0  =  R  U  - 1  ‘  '  . .  ‘  ”  . .  . '  1/  3 

_ U0  =  (  VU*VK1  )  D  I  *UKl _ 1  /  C 

M0VE(N)N0S.(-H0M(Vu,N/-1)TC  (V1.N/-1)  1  // 

_ mQvE(N)Nos.mum(Hq>n/-i)ti  xfi;L»N/'"l) _ _  _ _ _  i/b 

MOVE (N)NOS. mum (Uo/nz-dtc  <ci/n/-i)  .  -  -  i / v 


35 


f'0VHNJNOS.FHUM(U0,  N/-l)Tt  {U1.N/-1  ) 

......  16U 

161 

ibi 

V1  =  V0X  Ql*UQ*  K1»RoX  U 1 « U 0 

0010(0.8) 

Vi 

F OHH l 4-6) 1-12 ) 

Ibw4 

►  Ml 

FOHhilii-10)4.6) 

164 

►  F 

K)«M  (12-6)3-1)1-10) 

Id!? 

F  m2 

(•  oKMie-io  )b-io  > 

lod 

► 

►  ORM('fl-lOTb-4)<t-6)3-2)i-i)i2_ro)3-l)l-5)<*-6)" 

10  / 

LIST 

166 

LNl)  OQ  1  Q ( S I AR  T  ) 

lov 

BUUNDAHY  P I t 

.  yn 

.01  .  2t>  ,25  15.0 

0 

30000000  00  1UU0000U  01  0.10  i.oouooouoo 

5000Q 

RRU0 

REFERENCES 


1.  M.  J.  Romanelli,  "Introductory  Programming  for  ORDVAC  and  BRLESC," 
FORAST  (Formula  and  Assembly  Translator) ,  Ballistic  Research 
Laboratories  Report  No.  1209 ,  July  1963. 

2.  L.  W.  Campbell  and  G.  A.  Beck,  "The  FORAST  Programming  Language  for 
ORDVAC  and  BRLESC  (Revised),"  Ballistic  Research  Laboratories  Report 
No.  1273,  March  1965.  (Supersedes  Report  No.  1172). 


37 


DISTRIBUTION  LIST 


No.  of 

Copies  Organization 

20  Commander 

Defense  Documentation  Center 
ATTN :  TIPCR 
Cameron  Station 
Alexandria,  Virginia  22314 

1  Commanding  General 

U.S.  Army  Materiel  Command 
ATTN :  AMCRD-TE 

Washington,  D.C.  20315 

1  Commanding  General 

U.S.  Army  Missile  Command 
ATTN :  AMSMI-AML 

Redstone  Arsenal,  Alabama 
35809 

1  Commanding  Officer 

U.S.  Army  Engineer  Research 
&  Development  Laboratories 
ATTN:  STINFO  Div 

Fort  Belvoir,  Virginia  22060 

1  Commanding  Officer 

U.S.  Army  Picatinny  Arsenal 
ATTN :  SMUPA-VA6 

Dover,  New  Jersey  07801 

1  Commanding  General 

U.S.  Army  Weapons  Command 
Rock  Island,  Illinois  61202 

1  Commanding  Officer 

U.S.  Army  Rock  Island  Arsenal 
Rock  Island,  Illinois  61202 

1  Commanding  Officer 

U.S.  Army  Materials  Research 
Agency 

Watertown,  Massachusetts  02172 

1  Commanding  Officer 

U.S.  Army  Research  Office 
(Durham) 

Box  CM,  Duke  Station 
Durham,  North  Carolina  27706 


No.  of 

Copies  Organizat ion 
3  Commander 

U.S.  Naval  Air  Systems  Command 
Headquarters 
ATTN :  AIR-604 

Washington,  D.C.  20360 

1  Commander 

U.S.  Naval  Ordnance  Laboratory 
Silver  Spring,  Maryland  20910 

1  Superintendent 

U.S.  Naval  Postgraduate  School 
ATTN:  Tech  Rept  Sec 
Monterey,  California  93940 

1  Director 

U.S.  Naval  Research  Laboratory 
Washington,  D.C.  20390 

1  Commander 

U.S.  Naval  Weapons  Laboratory 
ATTN:  Tech  Lib,  MAL 

Dahlgren,  Virginia  22448 

1  RTD  (RTTM) 

Bolling  AFB,  D.C. 

20332 

1  AFML  (MAA) 

Wright-Patterson  AFB 
Ohio  45433 

1  ARL  (ARM,  Mr.  H.  Fettis) 
Wright-Patterson  AFB 
Ohio  45433 

1  Sandia  Corporation 

ATTN:  Org  1115,  Dr.  C.  Karnes 

P.0.  Box  5800 

Albuquerque,  New  Mexico  87115 

1  Sandia  Corporation 
Livermore  Laboratory 
ATTN:  Dr.  S.  Chiu 

P.0.  Box  969 

Livermore,  California  94550 


39 


DISTRIBUTION  LIST 


No.  of 

Copies  Organization 

2  General  Electric  Company 
Space  Sciences  Laboratory 
ATTN:  Dr.  F.  Wendt 

Dr.  J.  Heyda 
Dr.  T.  Riney 
P.  0.  Box  8555 
Philadelphia,  Pa.  19101 

3  Brown  University 
Division  of  Engineering 
ATTN:  Prof  P.  Symonds 

Prof  J.  Duffy 
Prof  D.  Drucker 
Providence,  Rhode  Island  02912  1 

3  Brown  University 

Division  of  Engineering 
ATTN:  Prof  H.  Kolsky 
Prof  R.  Clifton 
Prof  W.  Green 

Providence,  Rhode  Island  02912 

2  Cornell  University 
Department  of  Theoretical  & 

Applied  Mechanics 
ATTN:  Prof  D.  Robinson 

Prof  R.  Lance 
Ithaca,  New  York  14850 

3  The  University  of  Texas 
Department  of  Engineering 

Mechanics 

ATTN:  Prof  E.  Ripperger 

Prof  H.  Calvit 
Prof  C.  Yew 
Austin,  Texas  78758 

1  Professor  J.  Bell 

Department  of  Engineering 
Mechanics 

The  Johns  Hopkins  University 
34th  and  Charles  Street 
Baltimore,  Maryland  21218 


No.  of 

Copies  Organization 

1  Professor  H.  Bleich 

Institute  of  Flight  Structures 
Columbia  University 
New  York,  New  York  10027 

1  Professor  P.  C.  Chou 

Department  of  Aeronautical 
Engineering 

Drexel  Institute  of  Technology 
32nd  and  Chestnut  Streets 
Philadelphia,  Pa.  19104 

Professor  N.  Davids 
Department  of  Engineering 
Mechanics 

Pennsylvania  State  University 
319  Electrical  Engineering  East 
University  Park,  Pa.  16802 

1  Professor  W.  Goldsmith 
Department  of  Engineering 
Mechanics 

University  of  California 
Berkeley,  California  94704 

1  Professor  P.  D.  Lax 

Courant  Institute  of  Mathematical 
Sciences 

New  York  University 
251  Mercer  Street 
New  York,  New  York  10012 

1  Professor  E.  H.  Lee 

Department  of  Engineering 
Mechanics 

Stanford  University 
Palo  Alto,  California  9430H 

1  Professor  L.  Malvern 
College  of  Engineering 
Michigan  State  University 
East  Lansing,  Michigan  48823 


40 


DISTRIBUTION  LIST 


1  Professor  W.  Prager 

University  of  California 
8602  La  Jolla  Shores  Drive 
La  Jolla,  California  92037 


1  Professor  T.  Ting 

Department  of  Materials 
Engineering 
University  of  Chicago 
Chicago,  Illinois  60637 


1  Dr  .  T .  D .  Dudderar 

Bell  Telephone  Laboratory 

Mountain  Avenue 

Murray  Hill,  New  Jersey  07971 


1  Dr.  U.  Lindholm 

Department  of  Mechanical 
Sciences 

Southeast  Research  Institute 

8500  Culebra  Road 

San  Antonio,  Texas  78228 


1 


Dr.  R.  Sedney 

Research  Institute  for  Advanced 
Studies 

The  Martin  Company 
7212  Bellona  Avenue 
Ruxton,  Maryland  21212 


Aberdeen  Proving  Ground 


Ch,  Tech  Lib 

Air  Force  Ln  Ofc 
Marine  Corps  Ln  Ofc 
Navy  Ln  Ofc 
CDC  Ln  Ofc 


Ul 


Security  Classification 


DOCUMENT  CONTROL  DATA  -RAD 

(Security  c!aa«/fic«(ion  ot  tltlo,  body  ot  mbmtrmct  and  Indoxlng  annotation  muot  bo  mntorod  whon  tho  ororatt  roport  l»  eloomltlod 


1.  ORIGIN  A  TING  ACTIVITY  ( Cotpomt O  AtftflOfj  2 A.  REPORT  SECURITY  CLASSIFICATION 

ttoa  ,  T  ,  Unclassified 

U.S.  Army  Ballistic  Research  Laboratories  1 

Aberdeen  Proving  Ground,  Maryland 


3.  REPORT  TITLE 

THE  DYNAMIC  EXPANSION  OF  A  SPHERICAL  CAVITY  IN  AN  ELASTIC-PERFECTLY-PLASTIC  MATERIAL 


.  DESCRIPTIVE  NOTES  (Typo  ot  report  and  tnctuolvo  dm  tom) 


5.  AUTHOR(I)  (Firmt  nomo,  middlo  Initial,  loot  nmmo) 


nun. , 


«  REPORT  DATE 


February  1967 


l«.  CONTRACT  OR  GRANT  NO. 


7 A.  TOTAL  NO.  OP  PAGES  |7fr.  NO.  OP  PEPS 


•a.  ORIGINATOR’S  REPORT  NUMBER(S) 


b.  project  no. RDT&E  1P222901A201 


Report  No.  1357 


mb.  other  REPORT  NOISI  (Any  othor  numborm  that  moy  bo  ammlgpod 
thlm  roport) 


10.  DISTRIBUTION  STATEMENT 


Distribution  of  this  document  is  unlimited. 


12.  SPONSORING  MILITARY  ACTIVITY 

U.S.  Army  Materiel  Command 
Washington,  D.C. 


13.  ABSTRACT 


It  has  been  shown  that  a  finite-difference  numerical  technique  can  be  used  to  solve 
mixed  initial -and  boundary-value  problems  involving  high-speed  elastic-plastic  flow 
with  spherical  symmetry.  Numerical  solutions  for  the  dynamic  expansion  of  a 
spherical  cavity  under  a  constant  pressure  are  presented  to  demonstrate  the  nature 
and  capability  of  the  numerical  scheme.  The  solution  for  an  elastic  material  agrees 
closely  with  the  exact  solution.  The  solution  for  an  elastic-perf ectly-plastic  material 
has  confirmed  Green’s  prediction  concerning  the  motion  of  the  elastic-plastic  boundary. 
At  large  times,  the  asymptotic  solution  of  the  dynamic  problem  is  different  from  the 
quasi-static  solution.  This  result  indicates  that  the  quasi-static  approximation  may 
not  hold  in  dynamic  plasticity.  A  non-linear  dependence  of  the  plastic  solution  on 
the  boundary  condition  is  also  observed  in  the  results. 


REPLACES  DO  FORM  I47t,  I  J AN  M,  WHICH  IB 
OBSOLETE  POR  ARMY  USE. 


Unclassified 


Security  Classification 


